Field-Induced magnetic structures in Tb 2 Ti 2 07 at low temperatures: 
From spin-ice to spin-flip structures. 
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We studied the field-induced magnetic structures of the Tb2Ti207 pyrochlore by single-crystal 
neutron diffraction with a magnetic field applied along a [110] axis, focusing on the influence of a 
small misalignment. Both induced magnetic structures with k — and k — (0,0,1) propagation 
vectors are found to be sensitive to the misalignment, which controls the magnitude and orientation 
of the Tb moments involved in the f3 chains, with local [111] anisotropy axis perpendicular to the 
field. For k = 0, spin-ice-like structures are observed for a misalignment of a few degrees, whereas 
other structures, where the Tb-/3 moments flip by "melting" on the field axis, occur when the field 
is perfectly aligned. The field evolution of the k = structure is well reproduced by a molecular 
field model with anisotropic exchange. We give a complete symmetry analysis of the k — and 
k = (0,0, 1) magnetic structures, both being described by the basis functions of single irreducible 
representations. 
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I. INTRODUCTION 



Geometrically frustrated magnets such as spin liquids 
and spin ices bring out fundamental questions about their 
ground state and allow low-energy spin excitations to be 
observed. Historically 1 , a spin liquid is defined as a coop- 
erative paramagnet, where spins fluctuate down to zero 
temperature with spin-spin correlations rapidly decaying 
with increasing interatomic distances. From theory, it is 
now well established that the antiferromagnetic (AFM) 
interactions between isotropic (Heisenberg) spins yield a 
spin-liquid ground state for peculiar geometries, like the 
kagome, checkerboard, and pyrochlore lattices 1-3 . Quan- 
tum spin-liquid states have also been found in dimer 
models 4 , when the ground state is described by a su- 
perposition of spin singlets without symmetry breaking. 
The characteristic of spin liquids is the continuous de- 
generacy of the ground-state manifold, leading to macro- 
scopic entropy at T = OK. By contrast, spin ices refer to 
systems with discrete degeneracy issuing from topological 
constraints, like those which govern the two energetically 
equivalent positions of protons of water molecules in real 
ice. Hence, spin ices have the same residual entropy as 
water ice, and they have been found in kagome and py- 
rochlore lattices, where ferromagnetically coupled Ising 
spins are constrained to lie along their local anisotropy 
axes 5 ' 6 . Such constraints yield exotic excitations akin to 
magnetic monopoles 7 , which locally violate the ice rules 
by stochastic propagations along Dirac strings. 



Rare-earth titanate pyrochlores -R2T12O7, where the 
rare-earth ions R belong to a lattice of corner shar- 
ing tetrahedra, show well-known examples of spin-liquid 
(R = Tb) and spin-ice (R = Dy, Ho) behaviors. 
Tb2Ti207 is especially interesting since the origin of its 
spin-liquid ground state 8 remained mysterious for many 
years. Simple considerations would predict a transi- 
tion toward a nonfrustrated Neel order with Ising local 
anisotropy, in contrast to experiment. Indeed, Tb2Ti207 
shows dynamic short-range correlations, without onset of 
long-range magnetic order down to the lowest measured 
temperature of 50 mK. The energy terms responsible for 
this anomaly could be the peculiar crystal field scheme 
of the Tb 3+ ion, with a low-energy splitting between the 
ground and first excited doublets 9,10 , and an effective 
Tb-Tb first neighbor interaction at the verge between 
AFM and ferromagnetic. Recent theories suggest that 
crystal field fluctuations turn the effective Tb-Tb inter- 
action to ferromagnetic, leading to a fluctuating quantum 
spin-ice state 11 . Tb2Ti207 may therefore represent the 
"missing link" between spin liquid and spin ice, having a 
highly tunable ground state, easily destabilized by mag- 
netic field 12 , pressure 13 ' 14 , or chemical substitution 15 . 

In Tb2Ti207 magnetic long-range order can be in- 
duced by a magnetic field applied along the [110] direc- 
tion 12 ' 16 . The field-induced structures are of two types, 
with k = and k = (0,0,1) propagation vectors. The 
k = structure has non-zero magnetization, whereas 
the k — (0, 0, 1) one is antiferromagnetic. In our previ- 
ous work 16 ' 17 , we studied their local spin configurations 
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by neutron diffraction versus field and temperature. In 
the k = structure, we found that the local order in 
low fields and at low temperatures is akin to a spin ice, 
but with different values of the magnetic moments in a 
tetrahedron. The k = (0, 0, 1) structure was observed 
with a finite length scale, in a limited range of the phase 
diagram above 2 T and below 2 K, coexisting with the 
k = one. 

Recent measurements in Ho2Ti2C>7 spin ice have shown 
that the field-induced structures are very sensitive to the 
precise alignment of the field with respect to the crystal 
axes 18 . This effect was qualitatively observed by check- 
ing the field dependence of some magnetic Bragg peaks, 
and it was even used as a trick to tune the concentra- 
tion of magnetic monopoles or the dimer concentration in 
the Kasteleyn state 19-21 . It results from the huge uniax- 
ial anisotropy of the paramagnetic susceptibility in spin 
ices 22 . Considering the uniaxial anisotropy, when H is 
along the [110] axis, the field induced structures involve 
the so-called a and (3 chains 16 ' 23 along perpendicular di- 
rections. The a chains have their local anisotropy axis at 
36° from the field whereas the j3 chains have their easy 
axis perpendicular to H. In Ho2Ti 2 07, the k — (0,0, 1) 
structure was observed only when the field was perfectly 
aligned along the [110] crystal axis 18 ' 24 . 

Here, we have checked the influence of a field misalign- 
ment in Tb2Ti2C>7. By varying the field alignment with 
respect to a [110] axis in a systematic way, and refin- 
ing the magnetic structures in each case, we determine 
how the local spin arrangements are affected. The mis- 
alignment mainly influences the moments on the Tb-/3 
chains and the balance between the two structures. The 
k = (0, 0, 1) structure, which involves almost only the 
Tb-/3 moments, is favored when the alignment is opti- 
mal, and disappears when it exceeds 4.5°. This explains 
previous discrepancies concerning the transition line in 
the magnetic phase diagram of Tb2Ti20y. In the k = 
structure, the magnitude and orientation of the Tb-/3 mo- 
ments also strongly depend on the misalignment. The 
spin-ice orders observed in low fields are actually stable 
only when the field is slightly misaligned. An optimal 
alignment favors Tb-/3 moments opposite to the field, up 
to a compensation field of 1.5(5) T where they vanish, and 
above which they start to grow parallel to the field. All 
these features are quantitatively explained using a four 
sites molecular field model with anisotropic exchange. 

We performed a symmetry analysis, taking into ac- 
count the symmetry lowering with respect to the Fdim 
space group induced by a field applied along [110]. We 
find that both k = and k = (0, 0, 1) structures can 
be described by irreducible representations in the space 
group I4i/amd, the highest subgroup of Fd3m allowing 
the presence of a homogeneous magnetization M z along 
a twofold axis, invariant by symmetry. The symmetry 
analysis describes the magnetic structures with a limited 
number of parameters, including the case of a misoriented 
field. It provides a sound basis for the description of the 
pyrochlore lattice with [110] magnetic field by means of 



the a and /3 chains, which can be useful for real spin ices. 



II. EXPERIMENT 

A single crystal of Tb2Ti207 was grown from a sintered 
rod of the same nominal composition by the floating- 
zone technique, using a mirror furnace 16 . Neutron- 
diffraction studies were performed on the diffractomcter 
5C1 (A = 0.845 A) at the Orphee reactor of the Lab- 
oratoire Leon Brillouin, Saclay and on the CRG-CEA 
diffractometer D23 (A = 1.2815 A) at the Institut Laue- 
Langevin, Grenoble. We used unpolarized neutrons, col- 
lecting typically from 200 to 400 reflections for each data 
set. We collected data under external field applied close 
to a [110] direction. The best orientation we achieved cor- 
responds to a field misaligned from [110] by A<j) = 0.4° 
in the (001) plane to within an accuracy of ±0.3°. The 
field misalignment angle as well as the precision of the 
alignment were estimated based on the least-squares re- 
finement of the orientation matrix using 20 reflections 
measured several times. Then the field was progressively 
inclined by manual adjustment of the goniometer head 
from the [110] to the [110] direction thus always remain- 
ing in the (001) plane. The field misalignment angle A(j> 
was varied between —3° and +16°, the temperature T 
between 0.1K and 4K and the field H between 0T and 
12 T. The program FullProf (Ref. 25) was used to re- 
fine the crystal and magnetic structure parameters. 

Prior to the low-temperature studies with applied field 
the crystal was characterized in zero field. A total of 
289 reflections with sin6*/A < 0.55 A -1 were measured 
at 1.6 K and 44 unique reflections were obtained by av- 
eraging equivalents (i?i nt = 0.043) using the cubic space 
group Fdim. These were used to refine one oxygen po- 
sitional parameter, all the isotropic temperature factors, 
the scale and the extinction parameters (Rp — 0.037). 

III. RESULTS AND DISCUSSION 

A. Influence of the field misalignment on the phase 
diagram 

To determine the conditions of stability of the k = 
(0,0,1) structure, three Bragg reflections of different 
types were considered. The first reflection (113) has both 
nuclear and k = magnetic components. The other 
two reflections (002) and (112) are purely magnetic and 
belong to the k = and k = (0, 0, 1) structures, respec- 
tively. They appear only under an applied magnetic field. 
In Fig. 1, the intensities of these reflections are plotted 
at 1.6 K versus the magnetic field for different angles be- 
tween the field and the [110] direction. The intensity of 
the (113) peak is only weakly affected by the misalign- 
ment. In contrast, the intensities of the (002) and (112) 
Bragg peaks strongly vary in opposite ways. The (112) 
peak which has about the same intensity as the (002) 
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FIG. 1. (Color online) Field dependences of selected Bragg 
reflections of both magnetic structures with fc — and fc = 
(0,0,1) at 1.6 K and different field misalignment angles A<j) 
with respect to the [110] direction. Error bars are smaller 
than the symbol size if not given. The dotted lines are guide 
to the eyes. 



FIG. 3. (Color online) H-Tafm phase diagram at different 
field misalignment angles with respect to the [110] direction 
according to our present measurements (squares), Ref. 16 (cir- 
cles) and Ref. 12 (diamonds). The typical alignment error is 
less than 0.5°. The dotted lines are guide to the eyes. 



XlO 2 









■•».. 




1 


6 


(113) 






■■•o 


o 

fc = 


4 










o 




(002) 




□•■ 




fc = 


2 




*C3 [f]..- 











(112) 


P 


V 


fc = 
■■<>■ 


(0,0,1) 

o 



-15 -10 -5 5 10 
(deg) 



FIG. 2. (Color online) Intensities of selected reflections mea- 
sured at 1.6 K and 7T as a function of the field misalignment 
angle Ac/>. Error bars are smaller than the symbol size if not 
given. The dotted lines are guide to the eyes. 



peak for a well-aligned sample at 7T, is completely sup- 
pressed for a misalignment above 4.5°, whereas the in- 
tensity of the (002) peak increases by about a factor of 
2. Similar observations can be made from Fig. 2, where 
the intensities of the three Bragg peaks at 1.6 K and 7T 
are plotted versus the field misalignment angle A<j>. 

The suppression of the (112) peak and of all other 
peaks of the same family shows the vanishing of the 
fc = (0, 0, 1) structure. By measuring the temperature 
and field dependences of such peaks for a given field mis- 



alignment, we have drawn in Fig. 3 the transition lines 
in the H-Tafm phase diagram of Tb2Ti20y. For a mis- 
alignment of 1.5°, the fc = (0,0,1) structure persists up 
to the highest measured field, and its transition tempera- 
ture increases with the field, reaching 4 K at 12 T. Above 
4.5° misalignment, no fc = (0,0, 1) structure is observed. 
Comparison with previous results with other field mis- 
alignments from Refs. 12 and 16 shows that the onset of 
the fc = (0,0, 1) structure always occurs above a critical 
field of 2 T but that its stability range strongly depends 
on the misalignment. For a given field, the transition 
temperature decreases with increasing the misalignment. 
This explains the difference between the transition lines 
in H-Tafm phase diagram of Tb2Ti2C>7 previously re- 
ported in the literature 12 ' 16 . 

In the same way, the reentrant character of the fc 
: ( ) . 0, 1) structure quoted in Refs. 16 and 17, can now be 
attributed to the influence of the field misalignment. In 
these papers, the fc = (0, 0, 1) structure was found to be 
stabilized below 1 K in a limited field range, with an onset 
at 2 T and a maximum of the intensity of the AFM peaks 
and correlation length around 5T. Therefore, it was sug- 
gested that the fc = (0, 0, 1) structure would disappear 
at high fields. The present measurements performed at 
0.3 K up to very high fields show that this structure in- 
deed occurs at 2 T, but remains stable at least up to 12 T 
when the sample is well aligned. 

From the opposite variations in the (002) and (112) 
Bragg peaks, we also conclude that the balance between 
fc = and fc = (0, 0, 1) structures is affected by the 
field misalignment. Magnetic refinement and symmetry 
analysis, performed separately for each structure, yield a 
quantitative characterization of this interplay. 
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B. Symmetry analysis of the k = structures 

To analyze the field-induced magnetic structure of 
Tb 2 Ti 2 07, a systematic search was performed based on 
the theory of representations of space groups proposed by 
Bertaut 26 and Izyumov 27 . It allows one to consider all 
possible models of magnetic structures consistent with a 
given crystal structure of space group G. According to 
this method the magnetic structure can be expressed via 
the basis vectors of the irreducible representations of the 
group G. To calculate the irreducible representations the 
program Baslreps from the FullProf suite 25 was used. 

In zero field there are four irreducible representations 
for the space group Fd3m associated with the propa- 
gation vector k = 0. However, none of the zero-field 
magnetic structures allowed by the Fd3m symmetry are 
compatible with the experimental data obtained under 
field H || [110]. This suggests that the magnetic field 
applied along the twofold [110] axis breaks the Fd3m 
symmetry. Therefore, we searched for new solutions in 
the space group I4i/amd, the highest subgroup of Fd3m, 
for which a homogeneous magnetization component M z 
induced by the field applied along the two-fold axis is 
invariant. 

The crystal structure of i? 2 Ti 2 C>7 pyrochlore com- 
pounds in the Fdim space group is shown in Fig. 4 to- 
gether with the transformation from the Fd3m to the 
IAi/amd space group. We notice that although a dis- 
tortion of the cubic structure was inferred in Tb 2 Ti 2 07 
below 20 K from x-ray measurements 28 , it remains much 
too small to be observed with the resolution of neutron 
diffraction. The new unit-cell parameters in IA\/amd 

are a* = and c* = a, where a is the cubic unit-cell 
parameter of the FdZm space group. Hereafter, the sym- 
bol * will be used to denote the lattice parameters, coor- 
dinates and directions associated with the IAi/amd sym- 
metry. The positions of the four Tb sites (in IAi/amd) 
are defined in Table I. They are labeled according to the 
a (Tb-ai and Tb-a 2 ) or /3 (Tb-ft and Tb-/3 2 ) chains to 
which they belong. 

We start by considering that, in spite of the fact that 
the a and j3 chains can have different moments, they are 
both involved in the phase transition with k = 0. In this 
case it is generally accepted that their basis functions 
should belong to the same irreducible representation. 
The magnetic representation for the Tb 3+ ions occupying 
the 8d site in the I4i /amd space group can be written as 

r Tb = irf ) + ir^ + 2t£> + 2i^ 1) + 3r^ 2) , following 

the numbering scheme of Kovalev 29 , where the super- 
script is the order of the irreducible representation. The 
basis vectors describing the Tb moments which trans- 
form according to the irreducible representations r 1; r 3 , 
Tq and can be ruled out on the basis of magnetiza- 
tion measurements as none of them allows the existence 
of a net ferromagnetic moment along the field direction 
[010]*. In contrast, the two-dimensional representation 
Fg described in Table I allows the presence of both a 



TABLE I. Irreducible representation Tq of Tb2Ti207 (space 
group 1 4i /amd) associated with k — 0. Basis vectors pro- 
jected from a general vector M with the components M x , 
My, and M z at the Tb 8d sites. Tb-Qi: (x,y,z)*, Tb-a 2 : 
{-x + \,-V,z+ §)*, Tb-ft: {y+±,x+l,-z+ |)* and 
Th-/3 2 : (-y+ \,-x+ \,-z+ §)*. 
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ferromagnetic order along the [010]* direction (see the 
My component in Table I) and an antiferromagnetic or- 
der along the [001]* direction (see the M z component). 
Moreover, Tg also allows a net magnetization component 
along [100]* (see the M x component). As a result, the 
magnetic structure induced by a field misaligned from 
the [010]* direction but remaining in the (001)* plane 
can also be described by the same irreducible represen- 
tation Tg. 

We note that in the basis functions of Tg, the mag- 
netic moments in the Tb-a and j3 chains are completely 
decoupled and can be varied separately. To summarize, 
for the refinement of the k = structure one can use 
three fitting parameters U x , U y , and U z for the Tb-a 
chains and three other parameters V x ,V y , and V z for the 
Tb-/3 chains. Therefore, the symmetry analysis reduces 
the number of the refined parameters from twelve in an 
unconstrained refinement (three for each of the four Tb 
moments in a tetrahedron), to six only. 



C. Magnetic refinements of the k = structure 

To refine the k = structure of Tb 2 Ti 2 7 , from 200 
to 400 Bragg reflections were collected up to sin6>/A « 
0.55 A -1 . Such a collection was made in each case, 
namely, for given values of the magnetic field, temper- 
ature and field misalignment. The integrated intensities 
of these reflections were used to refine the components 
of the Tb 3+ magnetic moments using the symmetry con- 
straints described above. We note that with this reduced 
number of parameters, we obtain the same or even bet- 
ter quality of refinement as in previous unconstrained re- 
finements 16,17 . Nuclear structure parameters were taken 
from the low-temperature measurements in zero field. We 
found that, among all possible irreducible representations 
for k = 0, only T 9 (see Table I) provides a reliable de- 
scription of the neutron-diffraction data. 

Figure 5 shows the angles tj> and 9 of the Tb mag- 
netic moments of Tb 2 Ti 2 0y and typical moment config- 
urations for the k = magnetic structures are shown in 
Figs. 6 and 7. In all figures, quantities related to a and 
ft moments are shown by different colors, using the same 
conventions for clarity. The angles <fi and 9 refer to the 
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FIG. 4. (Color online) Left: crystal structure of Tb2Ti207. Only the Tb ions are shown for simplicity. Some tetrahedra are 
numbered for an easier comparison with Figs. 6, 7, 9 and 10. Right: transformation from the FdZm to IAi/amd space group; 
view along the [001] direction. 



[100] and [001] axes of the cubic unit cell, respectively, 
the field H is along the [110] axis. 

Figure 5 (left panel) shows the case of well-aligned 
fields (Acj) = 0.4° and 1.5°). As seen from the 9(H) de- 
pendence, when the field increases, the Tb-cv moments 
slowly rotate from their local [111] easy axis towards 
the field direction [110] lying in the (001) plane. Con- 
sequently, the antiferromagnetic component of the Tb-a 
moments along [001] becomes smaller. We note that at 
12 T the a moments still deviate, by about 12°, from the 
field direction. As for the Tb-/3 moments, we found that 
for a well-aligned field, they are collinear to the field. 
Surprisingly, they flip from a direction opposite to H (A 
state, Fig. 5) to a direction parallel to H (B state in 
Fig. 5), as shown by an abrupt change in the <fi angle. 
This flip of the Tb-/3 moments from an antiparallel to 
a parallel configuration occurs at a critical field H s { of 
1.5(5) T. Two typical spin structures below and above 
the critical field are drawn in Fig. 6. 

Figure 5 (right panel) shows the same quantities for 
fields misaligned (by A(p = 4.0° and 4.7°) toward the 
[100] direction. The comparison between the left and 
right panels clearly shows that the Tb-a moments are 
almost independent of the field misalignment. By con- 
trast, a significant change occurs for the /? chains. In the 
samples misaligned by 4.0° and 4.7° the field variations 
in the cp angle are much smoother showing that the Tb-/? 
moments rotate from their local [111] easy axis toward 
the (001) plane when the field increases. In such cases, 
a spin-ice-like configuration of the magnetic moments is 
stabilized in fields of about 1.5 T. This case corresponds 
to the spin-ice local order previously observed in Rcf. 16 
in the Tb2Ti 2 07 sample with a misalignment of 4° in the 
(001) plane, using unconstrained refinements. The spin- 



ice-likc configuration is illustrated in Fig. 7. We notice 
that the field region where the spin-ice local order is sta- 
bilized in the misaligned samples roughly coincides with 
the critical field H s f where the flip of the Tb-/3 moments 
occurs in the well-aligned samples. 

Another important information is given by the mo- 
ment values (Fig. 8). In the k — structure, the a mo- 
ments are practically independent of the field misalign- 
ment (Fig. 8, left panel) but the /3 moments strongly in- 
crease with the misalignment (Fig. 8, right panel). More- 
over, when the field is well aligned (Acf> — 0.4°), the /3 
moments show a minimum versus the field, decreasing 
to zero at -ff s f=1.5(5) T. This means that the spin flip 
occurs through a "melting" of the f3 moments. Namely, 
the j3 moments become paramagnetic at H s f because the 
applied field exactly compensates the exchange field. An- 
other example of spin "melting" is provided by the frus- 
trated XY model 30 . The vanishing of the /3 moments 
disappears together with the flip of the spins when the 
misalignment increases. One observes a small plateau for 
a misalignment of 4° , then a strong and monotonous in- 
crease when the field is strongly misaligned by 16°. In 
such a case, the magnitude and field dependence of the 
/? moments become comparable to those of the a ones, 
reflecting a gradual evolution from spin-ice-like to field- 
oriented structures with increasing field. 

The influence of a field misalignement with respect to 
[110] can be qualitatively understood by its decomposi- 
tion into a large [110] component along the a chains and a 
smaller perpendicular component, as suggested in Ref. 18 
for Ho2Ti207. In this case, the perpendicular component 
is along [110], namely, parallel to the j3 chains, and tilts 
the /3 moments. As to the a moments, a small misalign- 
ment should have little effect since the field component 
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FIG. 5. (Color online) Tb2Ti207 k — structure at 1.6 K. Field dependences of the Tb-a and ft moment angles for different field 
misalignments (A</>) with respect to the [110] direction. Error bars are smaller than the symbol size if not given. The angles cj> 
(top) and 9 (bottom) of a given moment /j, are defined at the center in spherical coordinates. The Q2 (/fe) moment is symmetric 
of Qi (/3i) with respect to the (001) plane. Left: well-oriented samples with A<f> = 0.4° and 1.5°. Right: misoriented samples 
with A(f> = 4° and 4.7°. The solid lines are calculations with the model described in Sec. HID. The following parameters were 
used for the calculations: Ax = — 0.42T//ib, An = — 0.06T//xb, A</> = 0° (left) and A<f> = 3° (right). The magnetic structures 
corresponding to the angles marked with "A", "-B" and "spin-ice-like" (vertical dotted lines) are shown in Figs. 6 and 7 and 
described in the text. 



along [110] decreases very slightly. A quantitative expla- 
nation of the field dependence of the k = structure is 
given by the model described below. 



D. Model calculation for the k = structure 

We implemented a model, based on the molecular 
field approximation, in order to reproduce the evolution 
with the field of the magnetic structure with k = 0. 
This model performs mean-field self-consistent calcula- 
tions and uses as ingredients the crystal field parame- 



ters of Tb 3 + in Tb 2 Ti 2 7 as in Refs. 10 and 22, and 
anisotropic two-ion exchange of the type: 

%ex = —J~\\ S\z S2Z — J±{Six S2X + Siy S2Y), (1) 

where ^ and ^7j_ are the components of the exchange 
tensor in the local frame with trigonal symmetry. 

We consider each ion in a tetrahedron to be exchange- 
coupled to its six nearest neighbors, resulting in an 
anisotropic molecular field tensor A such that the molec- 




FIG. 6. (Color online) Field-induced magnetic structure with 
k = of Tb 2 Ti 2 7 in the well-aligned sample (A(j> = 0.4°) 
below (A) and above (B) the flip of the Tb-/3 moments. Tetra- 
hedra are numbered for an easy comparison with Fig. 4. Mag- 
netic moment angles corresponding to structures denoted by 
A and B are shown in Fig. 5 and described in the text. 

ular field acting on this ion writes: 

1~ 6 

H mo \ = -\^2m k , (2) 

k=l 

where the sum runs over the six nearest neighbors. 

The relationship between the components of the A and 
J tensors is: 

X i =6J i (^ 1 ) i (3) 

where gj = 3/2 for Tb 3+ . 

A self-consistent treatment involving the four Tb mo- 
ments, each with its three components, is performed, the 
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FIG. 7. (Color online) Field-induced magnetic structure with 
k = of Tb 2 Ti 2 7 in the misaligned sample {A(j> = 4.7°). 
Tetrahedra are numbered for an easy comparison with Fig. 4. 
Such a spin-ice-like state is realized only for the misaligned 
fields of about 1.5 T. 



only parameters being thus the two components of the A 
tensor. This type of calculation is expected to describe 
only field induced magnetic structures having k = since 
only the moments on a single tetrahedron are considered. 
The dipole-dipole field is not included in the calculation. 

For perfect alignment of the field with [110], the model 
is able to reproduce quite well all the experimental data 
(see solid lines in Figs. 5 and 8) with an anisotropic an- 
tiferromagnetic A tensor. The spin-flip field H s f of the 
/? moments, in particular, is a rapidly increasing func- 
tion of the transverse component A_l and allows it to 
be determined with precision: Aj^ = —0.48(12) T///B, 
yielding H s f — 1.5(5) T. As to the longitudinal compo- 
nent An, we found that the isotropic exchange model 
with Aj^ = A|| does not match to the experiment. An 
agreement can be obtained only above a certain limit: 
A || > — 0.27T//ZB- In fact the calculated magnetic struc- 
tures show very weak dependence on An , but only above 
this value (— 0.27 T//xb) the calculated low-field magnetic 
structure match the experiment. 

The Aj_ value obtained here is about twice smaller than 
that (— 1T/hb) derived in Ref. 22 from the thermal de- 
pendence of the transverse susceptibility x_l(T). This 
can be due to the fact that the present model, involving 
four magnetic sub-lattices with six-neighbor exchange, 
is much more appropriate at low temperatures than the 
one-sublattice calculation performed in Ref. 22. Let us 
mention that calculation of x± (T) with the present model 
(not shown) yields very good agreement with the data 
using Aj^ = —0.42(8) T/^b, very close to the value de- 
rived above. This corresponds to a transverse exchange 
integral: J ± = -0.48 K. 

For the case of a misalignment of the magnetic field 




H (T) H (T) 

FIG. 8. (Color online) Tb2Ti207 fc — structure at 1.6 K. Field dependences of the Tb moment values for different field 
misalignments (Ac/>) with respect to the [110] direction. Error bars are smaller than the symbol size if not given. Left: Tb- 
a moments. Right: Tb-/3 moments. The solid lines are calculations with the model described in Sec. HID. The following 
parameters were used for the calculations: Ax = — 0.42T//ib, Aii = — 0.06 T/^b and Ac/> of 14°, 3° and 0° from top to bottom 
(right). The total j3 moment calculated as a vector sum of the fc = and fc = (0,0, 1) moments is also shown at the left. The 
dotted line is a guide to the eyes. 



with respect to [110], the model also reproduces quite 
well the data (see solid lines in Figs. 5 and 8) but the 
assessed misalignment angle in the (100) plane must be 
taken lower than the measured angle by about 0.5-1.5°. 
We notice that the results depend slightly on the sign of 
the misalignment, either tj> < 45° or <f> > 45°. 



E. Symmetry analysis and refinement of the 
k = (0, 0, 1) structures 

The symmetry analysis of the fc = (0, 0, 1) structure 
was performed in the space group I4i/amd as in the 
case of the fc = magnetic structure. In the same way, 
we find four two-dimensional irreducible representations 
(i — 1-4 according to Kovalev 29 notation) associated 

with k = (o, o, i), r^ b = i rf ) + 2 r' 2 (2) + 2 r' 3 (2) + 1 rf ) . 

Hereafter, the symbol ' denotes the parameters which 
correspond to the k = (0, 0, 1) structure. Among the 
representations, only T' 3 (see Table II) is in agreement 
with our experimental data. As seen from Table II four 
fitting parameters are needed to describe the T' 3 repre- 
sentation: two of them (U' and U' z ) describe the a chains 
and the remaining two (VJ. and V' z ) describe the (3 chains. 

The k = (0,0, 1) magnetic structure was investigated 
in the well-aligned samples under an applied field up to 
12 T. In the case of Acj) = 0.4°, a total of 196 pure mag- 
netic reflections with sin e/x < 0.55A- 1 were measured 
at 1.6 K and 7T. They were used to refine four Tb mag- 
netic moment components (Rf = 0.046) taking into ac- 
count the symmetry constraints of T' 3 . The magnetic 
structure resulting from the refinement is illustrated in 
Fig. 9. We find that the Tb-a moments are very small, 



TABLE II. Irreducible representation V' 3 of Tb2Ti2C>7 (space 
group IA\/amd) associated with fc = (0,0,1). Basis vectors 
projected from a general vector M with the components M x , 
My, and M z at the Tb 8d sites. The Tb positions are defined 
in Table I. 



Irr. Rep. 


Atom 


M x 


My 


M z 


r' 3 


Tb-ai 





U'y 


u' z 




Tb-a 2 







-u'. 




Tb-/?! 


v; : 





-v; 




Tb-/3 2 







K 



typically 0.2-0.6 ^b, as already quoted in Ref. 16, so that 
the magnetic order mostly concerns the Tb-/3 moments. 
The k = (0, 0, 1) propagation vector yields an AFM cou- 
pling of the P chains. Within one chain, the /3 moments 
order antiferromagnetically along the chains. We notice 
that such AFM order along the /3 chains contrasts with 
the k = X magnetic structure observed in rIo 2 Ti 2 07 
spin ice 24 . In Ho2Ti207, the (3 chains also couple anti- 
ferromagnetically but the moments within one chain are 
ferromagnetically coupled along the chain. 

F. Resulting magnetic structure of Tb2Ti2C>7 

First we recall that the fc = (0,0, 1) structure is stabi- 
lized only at low temperatures, in high and well-aligned 
fields (typically above 2 T and for a misalignment below 
4.5°). In all other cases, only the fc = structure is 
observed. In the (H, T, Ac/)) range where the two struc- 
tures coexist, one can easily reconstruct the full mag- 
netic structure of Tb-2Ti207 from the refined fc = and 
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k = (0,0,1) 
7T, 0.4 deg 



FIG. 9. (Color online) Field-induced magnetic structure with 
fc = (0,0, 1) of Tb 2 Ti 2 7 in the well-aligned sample (Acj> = 
0.4°). Tetrahedra are numbered for an easy comparison with 
Fig. 4. 




k = + 

k = (0,0,1) 
7T, 0.4 dee 



FIG. 10. (Color online) Field-induced magnetic structure of 
Tb2Ti207 in the well-aligned sample (A(j> = 0.4°). A super- 
position of fe = and fc = (0, 0, 1) is shown. Tetrahedra are 
numbered for an easier comparison with Fig. 4. 



fe = (0, 0, 1) structures, as a coherent superposition of the 
magnetic modes involved in Tg and T' 3 representations. 
This resulting magnetic structure of Tb2Ti207 at 7T is 
illustrated in Fig. 10. 

The field dependences of the Tb magnetic moments re- 
sulting from this superposition are shown in Fig. 8. The a 
moments which are involved in the fc = structure only 
increase with the field and reach 7/Zb at 12 T, whatever 
the field alignment. By contrast, the (3 moments, which 
share between the two structures, strongly depend on 
it. In well-orientated samples at high fields, a /3-moment 
component of about 4 hb belongs to the fc = struc- 



a. 



2 - 



H (T) 

FIG. 11. (Color online) Magnetization M of Tb 2 Ti 2 7 mea- 
sured at 5K versus the magnetic field H applied along the 
[110] axis according to Ref. 31 (circles). The solid line is 
a calculation with the model described in Sec. HID using 
the molecular field tensor deduced from the present measure- 
ments. 



ture and another component of about 4 fj,B goes to the 
fc = (0, 0, 1) structure. However, the total moment cal- 
culated as the vector sum of the fc = and fc = (0, 0, 1) 
moments (Fig. 8, left panel) reaches only about 5.5 Hb, 
thus remaining smaller than the a moment. When the 
misalignment is above 4.5° and the fc = (0, 0, 1) struc- 
ture is suppressed, the j3 moments in the fc = struc- 
ture reach the value of the a moments of about 7 fiB in 
high fields. We notice that this 7 hb value is above the 
value of 5.1 /zb found at H = 0T for the fluctuating Tb 
moments 10 . This is explained by the mixing of the low- 
energy crystal field levels by the applied field. 

From our measurements and the mean field calculation 
described above, we can also calculate the magnetization 
and compare it to its determination by bulk magnetic 
measurements 31 . As shown in Fig. 11, both quantities 
are in perfect agreement, which is an important check of 
the robustness and consistency of our analysis and which 
ensures that all the experimental corrections have been 
done properly. One also notices that misorientation ef- 
fects become negligible here since the bulk magnetization 
is almost independent of the crystal orientation 31 . 



IV. CONCLUSION 

To summarize, we studied the field-induced magnetic 
structures in Tb2Ti2C>7 under magnetic field close to the 
[110] axis, by varying the field alignment in a systematic 
way. We observed different types of non-collinear struc- 
tures which are extremely sensitive to the field alignment. 
This sensitivity mostly concerns the moments on the /3 
chains whose local anisotropy axis is perpendicular to the 
field. In low fields (H < 2 T), when only the fc = struc- 
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ture is stabilized, spin-ice local orders are favored for a 
misalignment of a few degrees. Local orders where the 
Tb-/3 moments flip along the field direction and melt at 
a compensation field, occur when the sample is very well 
aligned (within 0.4°). 

We note that the mechanism in the flip of the Tb-/3 
moments in Tb2Ti2C>7 is considerably different from the 
spin-flop and spin-flip transitions in conventional antifcr- 
romagnets such as MnF2 and FeCl2, respectively 32 ' 33 . In 
these materials, the moment reorientation occurs as a re- 
sult of the competition between anisotropy, exchange and 
Zeeman energy, it is accompanied by a symmetry break- 
ing and it is characterized by a moment conservation. By 
contrast, no symmetry breaking occurs in Tb2Ti207 at 
the critical field H s { since both magnetic structures above 
and below H s f belong to the same irreducible representa- 
tion. Moreover, the magnetic moment of the Tb-/3 chains 
gradually decreases until the applied field exactly com- 
pensates the exchange field for H = H s f. At this point, 
the {3 moments become paramagnetic, which opens the 
possibility to observe new kinds of field induced excita- 
tions in Tb2Ti207. The parent compound E^T^Oy re- 
cently studied offers an intermediate case, with a decrease 
of the Er moment values at the spin-flip transition 34 . 

The field behavior of the k — structure is well ex- 
plained by a molecular field model with anisotropic ex- 
change. In high fields (H > 2T), the antiferromagnetic 
structure with k — (0, 0, 1) propagation vector, which in- 
volves mostly the j3 moments, is favored in well-aligned 



samples and rapidly disappears in the misaligned sam- 
ples. In this structure, the ft moments are antiferro- 
magnetically coupled along the chains, in contrast with 
Ho2Ti2C>7 spin ice. The high sensitivity of the magnetic 
orders to the field alignment results from the combination 
of the crystal-field anisotropy and pyrochlore geometry, 
which tunes the coupling of the (3 chains to the applied 
field. 

We performed a full symmetry analysis of the field in- 
duced magnetic structures in the pyrochlore lattice. The 
field evolution of the magnetic structure in Tb2Ti 2 07 is 
successfully described using single irreducible represen- 
tations of the IAi/amd space group, both for the k = 
and k = (0,0, 1) structures. The main advantage of the 
symmetry analysis is a reduction of the fitted param- 
eters needed in the refinement, six parameters for the 
k = structure and only four for the k = (0, 0, 1) one, 
instead of twelve in the unconstrained refinement. Simi- 
lar approaches may be used to analyze the field-induced 
magnetic structures of other pyrochlores. 

Note added: Recently, we were made aware of the 
study from another group 35 on the same subject. 
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